Selection leads to remarkable variability in the outcomes of hybridisation across replicate hybrid zones.

Hybrid zones have been viewed as an opportunity to see speciation in action. When hybrid zones are replicated, it is assumed that if the same genetic incompatibilities are maintaining reproductive isolation across all instances of secondary contact, those incompatibilities should be identifiable by consistent patterns in the genome. In contrast, changes in allele frequencies due to genetic drift should be idiosyncratic for each hybrid zone. To test this assumption, we simulated 20 replicates of each of 12 hybrid zone scenarios with varied genetic incompatibilities, rates of migration, selection and different starting population size ratios of parental species. We found remarkable variability in the outcomes of hybridisation in replicate hybrid zones, particularly with Bateson-Dobzhansky-Muller incompatibilities and strong selection. We found substantial differences among replicates in the overall genomic composition of individuals, including admixture proportions, inter-specific ancestry complement and number of ancestry junctions. Additionally, we found substantial variation in genomic clines among replicates at focal loci, regardless of locus-specific selection. We conclude that processes other than selection are responsible for some consistent outcomes of hybridisation, whereas selection on incompatibilities can lead to genomically widespread and highly variable outcomes. We highlight the challenge of mapping between pattern and process in hybrid zones and call attention to how selection against incompatibilities will commonly lead to variable outcomes. We hope that this study informs future research on replicate hybrid zones and encourages further development of statistical techniques, theoretical models and exploration of additional axes of variation to understand reproductive isolation.

Temperature impacts the bovine ex vivo immune response towards Mycoplasmopsis bovis.

Although cattle are the mammalian species with most global biomass associated with a huge impact on our planet, their immune system remains poorly understood. Notably, the bovine immune system has peculiarities such as an overrepresentation of γδ T cells that requires particular attention, specifically in an infectious context. In line of 3R principles, we developed an ex vivo platform to dissect host-pathogen interactions. The experimental design was based on two independent complementary readouts: firstly, a novel 12-14 color multiparameter flow cytometry assay measuring maturation (modulation of cell surface marker expression) and activation (intracellular cytokine detection) of monocytes, conventional and plasmacytoid dendritic cells, natural killer cells, γδ T cells, B and T cells; secondly, a multiplex immunoassay monitoring bovine chemokine and cytokine secretion levels. The experiments were conducted on fresh primary bovine blood cells exposed to Mycoplasmopsis bovis (M. bovis), a major bovine respiratory pathogen. Besides reaffirming the tight cooperation of the different primary blood cells, we also identified novel key players such as strong IFN-γ secreting NK cells, whose role was so far largely overlooked. Additionally, we compared the host-pathogen interactions at different temperatures, including commonly used 37 °C, ruminant body temperature (38-38.5 °C) and fever (≥ 39.5 °C). Strikingly, working under ruminant physiological temperature influenced the capacity of most immune cell subsets to respond to M. bovis compared to 37 °C. Under fever-like temperature conditions the immune response was impaired compared to physiological temperature. Our experimental approach, phenotypically delineating the bovine immune system provided a thorough vision of the immune response towards M. bovis and the influence of temperature towards that immune response.

Increased heterozygosity in low-pass sequencing data allows identification of blood chimeras in cattle.

In about 90% of multiple pregnancies in cattle, shared blood circulation between fetuses leads to genetic chimerism in peripheral blood and can reduce reproductive performance in heterosexual co-twins. However, the early detection of heterosexual chimeras requires specialized tests. Here, we used low-pass sequencing data with a median coverage of 0.64× generated from blood samples of 322 F1 crosses between beef and dairy cattle and identified 20 putative blood chimeras through increased levels of genome-wide heterozygosity. In contrast, for 77 samples with routine SNP microarray data generated from hair bulbs of the same F1s, we found no evidence of chimerism, simultaneously observing high levels of genotype discordance with sequencing data. Fifteen out of 18 reported twins showed signs of blood chimerism, in line with previous reports, whereas the presence of five alleged singletons with strong signs of chimerism suggests that the in-utero death rate of co-twins is at the upper limit of former estimates. Together, our results show that low-pass sequencing data allow reliable screening for blood chimeras. They further affirm that blood is not recommended as a source of DNA for the detection of germline variants.

Local Adaptation and the Evolution of Genome Architecture in Threespine Stickleback.

Theory predicts that local adaptation should favor the evolution of a concentrated genetic architecture, where the alleles driving adaptive divergence are tightly clustered on chromosomes. Adaptation to marine versus freshwater environments in threespine stickleback has resulted in an architecture that seems consistent with this prediction: divergence among populations is mainly driven by a few genomic regions harboring multiple quantitative trait loci for environmentally adapted traits, as well as candidate genes with well-established phenotypic effects. One theory for the evolution of these "genomic islands" is that rearrangements remodel the genome to bring causal loci into tight proximity, but this has not been studied explicitly. We tested this theory using synteny analysis to identify micro- and macro-rearrangements in the stickleback genome and assess their potential involvement in the evolution of genomic islands. To identify rearrangements, we conducted a de novo assembly of the closely related tubesnout (Aulorhyncus flavidus) genome and compared this to the genomes of threespine stickleback and two other closely related species. We found that small rearrangements, within-chromosome duplications, and lineage-specific genes (LSGs) were enriched around genomic islands, and that all three chromosomes harboring large genomic islands have experienced macro-rearrangements. We also found that duplicates and micro-rearrangements are 9.9× and 2.9× more likely to involve genes differentially expressed between marine and freshwater genotypes. While not conclusive, these results are consistent with the explanation that strong divergent selection on candidate genes drove the recruitment of rearrangements to yield clusters of locally adaptive loci.

Comparing genome scans among species of the stickleback order reveals three different patterns of genetic diversity.

Comparing genome scans among species is a powerful approach for investigating the patterns left by evolutionary processes. In particular, this offers a way to detect candidate genes that drive convergent evolution. We compared genome scan results to investigate if patterns of genetic diversity and divergence are shared among divergent species within the stickleback order (Gasterosteiformes): the threespine stickleback (Gasterosteus aculeatus), ninespine stickleback (Pungitius pungitus), and tubesnout (Aulorhynchus flavidus). Populations were sampled from the southern and northern edges of each species' range, to identify patterns associated with latitudinal changes in genetic diversity. Weak correlations in genetic diversity (F ST and expected heterozygosity) and three different patterns in the genomic landscape were found among these species. Additionally, no candidate genes for convergent evolution were detected. This is a counterexample to the growing number of studies that have shown overlapping genetic patterns, demonstrating that genome scan comparisons can be noisy due to the effects of several interacting evolutionary forces.

Model-based genotype and ancestry estimation for potential hybrids with mixed-ploidy.

Non-random mating among individuals can lead to spatial clustering of genetically similar individuals and population stratification. This deviation from panmixia is commonly observed in natural populations. Consequently, individuals can have parentage in single populations or involving hybridization between differentiated populations. Accounting for this mixture and structure is important when mapping the genetics of traits and learning about the formative evolutionary processes that shape genetic variation among individuals and populations. Stratified genetic relatedness among individuals is commonly quantified using estimates of ancestry that are derived from a statistical model. Development of these models for polyploid and mixed-ploidy individuals and populations has lagged behind those for diploids. Here, we extend and test a hierarchical Bayesian model, called entropy, which can use low-depth sequence data to estimate genotype and ancestry parameters in autopolyploid and mixed-ploidy individuals (including sex chromosomes and autosomes within individuals). Our analysis of simulated data illustrated the trade-off between sequencing depth and genome coverage and found lower error associated with low-depth sequencing across a larger fraction of the genome than with high-depth sequencing across a smaller fraction of the genome. The model has high accuracy and sensitivity as verified with simulated data and through analysis of admixture among populations of diploid and tetraploid Arabidopsis arenosa.

Large-scale mutation in the evolution of a gene complex for cryptic coloration.

The types of mutations affecting adaptation in the wild are only beginning to be understood. In particular, whether structural changes shape adaptation by suppressing recombination or by creating new mutations is unresolved. Here, we show that multiple linked but recombining loci underlie cryptic color morphs of Timema chumash stick insects. In a related species, these loci are found in a region of suppressed recombination, forming a supergene. However, in seven species of Timema, we found that a megabase-size "supermutation" has deleted color loci in green morphs. Moreover, we found that balancing selection likely contributes more to maintaining this mutation than does introgression. Our results show how suppressed recombination and large-scale mutation can help to package gene complexes into discrete units of diversity such as morphs, ecotypes, or species.

Long-term balancing selection on chromosomal variants associated with crypsis in a stick insect.

How polymorphisms are maintained within populations over long periods of time remains debated, because genetic drift and various forms of selection are expected to reduce variation. Here, we study the genetic architecture and maintenance of phenotypic morphs that confer crypsis in Timema cristinae stick insects, combining phenotypic information and genotyping-by-sequencing data from 1,360 samples across 21 populations. We find two highly divergent chromosomal variants that span megabases of sequence and are associated with colour polymorphism. We show that these variants exhibit strongly reduced effective recombination, are geographically widespread and probably diverged millions of generations ago. We detect heterokaryotype excess and signs of balancing selection acting on these variants through the species' history. A third chromosomal variant in the same genomic region likely evolved more recently from one of the two colour variants and is associated with dorsal pattern polymorphism. Our results suggest that large-scale genetic variation associated with crypsis has been maintained for long periods of time by potentially complex processes of balancing selection.

Transitions between phases of genomic differentiation during stick-insect speciation.

Speciation can involve a transition from a few genetic loci that are resistant to gene flow to genome-wide differentiation. However, only limited data exist concerning this transition and the factors promoting it. Here, we study phases of speciation using data from >100 populations of 11 species of Timema stick insects. Consistent with early phases of genic speciation, adaptive colour-pattern loci reside in localized genetic regions of accentuated differentiation between populations experiencing gene flow. Transitions to genome-wide differentiation are also observed with gene flow, in association with differentiation in polygenic chemical traits affecting mate choice. Thus, intermediate phases of speciation are associated with genome-wide differentiation and mate choice, but not growth of a few genomic islands. We also find a gap in genomic differentiation between sympatric taxa that still exchange genes and those that do not, highlighting the association between differentiation and complete reproductive isolation. Our results suggest that substantial progress towards speciation may involve the alignment of multi-faceted aspects of differentiation.

Causes and consequences of large clonal assemblies in a poplar hybrid zone.

Asexual reproduction is a common and fundamental mode of reproduction in plants. Although persistence in adverse conditions underlies most known cases of clonal dominance, proximal genetic drivers remain unclear, in particular for populations dominated by a few large clones. In this study, we studied a clonal population of the riparian tree Populus alba in the Douro river basin (northwestern Iberian Peninsula) where it hybridizes with Populus tremula, a species that grows in highly contrasted ecological conditions. We used 73 nuclear microsatellites to test whether genomic background (species ancestry) is a relevant cause of clonal success, and to assess the evolutionary consequences of clonal dominance by a few genets. Additional genotyping-by-sequencing data were produced to estimate the age of the largest clones. We found that a few ancient (over a few thousand years old) and widespread genets dominate the population, both in terms of clone size and number of sexual offspring produced. Interestingly, large clones possessed two genomic regions introgressed from P. tremula, which may have favoured their spread under stressful environmental conditions. At the population level, the spread of large genets was accompanied by an overall ancient (>0.1 Myr) but soft decline of effective population size. Despite this decrease, and the high clonality and dominance of sexual reproduction by large clones, the Douro hybrid zone still displays considerable genetic diversity and low inbreeding. This suggests that even in extreme cases as in the Douro, asexual and sexual dominance of a few large, geographically extended individuals does not threaten population survival.

The genetic architecture of hybrid incompatibilities and their effect on barriers to introgression in secondary contact.

Genetic incompatibilities are an important component of reproductive isolation. Although theoretical studies have addressed their evolution, little is known about their maintenance when challenged by potentially high migration rates in secondary contact. Although theory predicts that recombination can erode barriers, many empirical systems have been found to retain species-specific differences despite substantial gene flow. By simulating whole genomes in individuals of hybridizing species, we find that the genetic architecture of two contrasting models of epistatic hybrid incompatibilities and the context of hybridization can substantially affect species integrity and genomic heterogeneity. In line with theory, our results show that intergenomic incompatibilities break down rapidly by recombination, but can maintain genome-wide differentiation under very limited conditions. By contrast, intragenomic interactions that arise from genetic pathways can maintain species-specific differences even with high migration rates and gene flow, whereas introgression at large parts of the genome can simultaneously remain extensive, consistent with empirical observations. We discuss the importance of intragenomic interactions in speciation and consider how this form of epistatic fitness variation is implicated and supported by other theoretical and empirical studies. We further address the relevance of replicates and knowledge of context when investigating the genomics of speciation.

Unexpected ancestry of Populus seedlings from a hybrid zone implies a large role for postzygotic selection in the maintenance of species.

In the context of potential interspecific gene flow, the integrity of species will be maintained by reproductive barriers that reduce genetic exchange, including traits associated with prezygotic isolation or poor performance of hybrids. Hybrid zones can be used to study the importance of different reproductive barriers, particularly when both parental species and hybrids occur in close spatial proximity. We investigated the importance of barriers to gene flow that act early vs. late in the life cycle of European Populus by quantifying the prevalence of homospecific and hybrid matings within a mosaic hybrid zone. We obtained genotypic data for 11 976 loci from progeny and their maternal parents and constructed a Bayesian model to estimate individual admixture proportions and hybrid classes for sampled trees and for the unsampled pollen parent. Matings that included one or two hybrid parents were common, resulting in admixture proportions of progeny that spanned the whole range of potential ancestries between the two parental species. This result contrasts strongly with the distribution of admixture proportions in adult trees, where intermediate hybrids and each of the parental species are separated into three discrete ancestry clusters. The existence of the full range of hybrids in seedlings is consistent with weak reproductive isolation early in the life cycle of Populus. Instead, a considerable amount of selection must take place between the seedling stage and maturity to remove many hybrid seedlings. Our results highlight that high hybridization rates and appreciable hybrid fitness do not necessarily conflict with the maintenance of species integrity.

Genomic scan for single nucleotide polymorphisms reveals patterns of divergence and gene flow between ecologically divergent species.

Recent advances in population genomics have triggered great interest in the genomic landscape of divergence in taxa with 'porous' species boundaries. One important obstable of previous studies of this topic was the low genomic coverage achieved. This issue can now be overcome by the use of 'next generation' or short-read DNA-sequencing approaches capable of assaying many thousands of single nucleotide polymorphisms (SNPs) in divergent species. We have scanned the 'porous' genomes of Populus alba and Populus tremula, two ecologically divergent hybridizing forest trees, using >38,000 SNPs assayed by restriction site associated DNA (RAD) sequencing. Windowed analyses indicate great variation in genetic divergence (e.g. the proportion of fixed SNPs) between species, and these results are unlikely to be strongly biased by genomic features of the Populus trichocarpa reference genome used for SNP calling. Divergence estimates were significantly autocorrelated (P < 0.01; Moran's I up to 0.6) along 11 of 19 chromosomes. Many of these autocorrelations involved low divergence blocks, thus suggesting that allele sharing was caused by recurrent gene flow rather than shared ancestral polymorphism. A conspicuous low divergence block of three megabases was detected on chromosome XIX, recently put forward as an incipient sex chromosome in Populus, and was largely congruent with introgression of mapped microsatellites in two natural hybrid zones (N > 400). Our results help explain the origin of the 'genomic mosaic' seen in these taxa with 'porous' genomes and suggest rampant introgression or extensive among-species conservation of an incipient plant sex chromosome. RAD sequencing holds great promise for detecting patterns of divergence and gene flow in highly divergent hybridizing species.

Recombinant hybrids retain heterozygosity at many loci: new insights into the genomics of reproductive isolation in Populus.

The maintenance of species barriers in the face of gene flow is often thought to result from strong selection against intermediate genotypes, thereby preserving genetic differentiation. Most speciation genomic studies thus aim to identify exceptionally divergent loci between populations, but divergence will be affected by many processes other than reproductive isolation (RI) and speciation. Through genomic studies of recombinant hybrids sampled in the wild, genetic variation associated with RI can be observed in situ, because selection against incompatible genotypes will leave detectable patterns of variation in the hybrid genomes. To better understand the mechanisms directly involved in RI, we investigated three natural 'replicate' hybrid zones between two divergent Populus species via locus-specific patterns of ancestry across recombinant hybrid genomes. As expected, genomic patterns in hybrids and their parental species were consistent with the presence of underdominant selection at several genomic regions. Surprisingly, many loci displayed greatly increased between-species heterozygosity in recombinant hybrids despite striking genetic differentiation between the parental genomes, the opposite of what would be expected with selection against intermediate genotypes. Only a limited, reproducible set of genotypic combinations was present in hybrid genomes across localities. In the absence of clearly delimited 'hybrid habitats', our results suggest that complex epistatic interactions within genomes play an important role in advanced stages of RI between these ecologically divergent forest trees. This calls for more genomic studies that test for unusual patterns of genomic ancestry in hybridizing species.

Ever-young sex chromosomes in European tree frogs.

Non-recombining sex chromosomes are expected to undergo evolutionary decay, ending up genetically degenerated, as has happened in birds and mammals. Why are then sex chromosomes so often homomorphic in cold-blooded vertebrates? One possible explanation is a high rate of turnover events, replacing master sex-determining genes by new ones on other chromosomes. An alternative is that X-Y similarity is maintained by occasional recombination events, occurring in sex-reversed XY females. Based on mitochondrial and nuclear gene sequences, we estimated the divergence times between European tree frogs (Hyla arborea, H. intermedia, and H. molleri) to the upper Miocene, about 5.4-7.1 million years ago. Sibship analyses of microsatellite polymorphisms revealed that all three species have the same pair of sex chromosomes, with complete absence of X-Y recombination in males. Despite this, sequences of sex-linked loci show no divergence between the X and Y chromosomes. In the phylogeny, the X and Y alleles cluster according to species, not in groups of gametologs. We conclude that sex-chromosome homomorphy in these tree frogs does not result from a recent turnover but is maintained over evolutionary timescales by occasional X-Y recombination. Seemingly young sex chromosomes may thus carry old-established sex-determining genes, a result at odds with the view that sex chromosomes necessarily decay until they are replaced. This raises intriguing perspectives regarding the evolutionary dynamics of sexually antagonistic genes and the mechanisms that control X-Y recombination.